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An  analytic  solution  of  a  thermocouple  has  been  developed  in  order  to  gain  a  deeper  understanding  of  the 
physics  of  a  real  device.  The  model  is  established  for  both  rectangular  and  cylindrical  couples  and  is  made 
to  account  for  thermal  resistance  of  the  hot  and  cold  shoes  and  lateral  heat  transfer.  A  set  of  dimension¬ 
less  parameters  have  been  developed  to  determine  couple  behavior  and  serve  as  simplifying  justifica¬ 
tions.  New  dimensionless  parameters.  Device  Design  Factor  and  Fin  Factor,  are  introduced  to  account 
for  the  thermal  resistance  and  lateral  heat  transfer,  respectively.  Design  guidelines  on  couple  length 
and  cross-sectional  area  have  been  established  to  account  for  conditions  encountered  by  a  realistic  cou¬ 
ple.  As  a  result  of  thermal  resistances  a  lower  limit  on  the  length  of  the  couple  can  be  established.  In  the 
case  of  a  lateral  heat  transfer  couple  the  efficiency  is  found  to  depend  upon  cross-sectional  area  of  the  leg 
in  such  a  fashion  as  to  suggest  the  need  to  design  large  area  couples.  The  classic  thermoelectric  solution 
neglects  the  effect  of  thermal  resistance  and  lateral  heat  transfer,  leading  to  an  over  estimated  conversion 
efficiency.  The  work  presented  provides  a  path  to  incorporate  these  neglected  factors  and  offers  a  simpli¬ 
fied  estimation  for  couple  performance  based  on  analytic  solutions  of  governing  equations. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Methodology 

The  foundation  of  the  thermoelectric  device  is  the  fundamental 
thermoelectric  couple,  composed  of  two  conducting  legs  classically 
configured  electrically  in  series  and  thermally  in  parallel.  The  cou¬ 
pled  transport  of  thermal  and  electrical  energy  through  the  legs 
serves  to  operate  as  a  heat  pump  or  heat  engine  depending  on 
the  driving  mechanism.  Therefore  it  is  of  critical  interest  to  have 
a  clear  understanding  of  the  physics  of  a  thermoelectric  couple  in 
order  to  design  both  materials  and  devices  capable  of  providing 
the  desired  thermoelectric  output  in  an  efficient  and  optimized 
manner. 
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The  governing  equations  of  a  couple  can  be  derived  from  apply¬ 
ing  the  generalized  thermoelectric  form  of  Ohm’s  and  Fourier’s 
Laws  to  a  control  volume,  obtainable  from  the  conjugate  flux  force 
relations  of  Onsager  [1,2],  A  similar  approach  can  be  found  in  the 
classic  conduction  equation  in  many  heat  transfer  texts,  such  as 
Arpaci  [3],  Ozisik  [4],  and  Carslaw  and  Jaeger  [5],  The  couples  were 
investigated  as  steady  state  one  dimensional  domains,  neglecting 
temperature  and  electrical  gradients  in  directions  other  than  the 
primary  leg  length,  and  with  constant  material  properties.  Steady 
state  is  justified  by  common  application  of  devices,  and  for  the  pur¬ 
pose  of  obtaining  concise  solutions.  Investigation  into  transient 
behavior  has  been  performed  by  Meng  et  al.  [6],  Alata  et  al.  [7], 
and  Montecucco  et  al.  [8]  among  others.  The  one  dimensional 
assumption  can  be  justified  by  investigating  Biot  numbers  for  the 
primary  leg  direction  versus  the  lateral  direction  of  typical  thermo¬ 
electric  legs.  The  constant  material  properties  simplification  was 
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used  to  form  a  linear  model  so  that  meaningful  solutions  could  be 
developed.  The  effect  of  variable  material  properties  is  out  of  the 
scope  of  this  work.  Some  of  the  effects  of  variable  material  proper¬ 
ties  may  be  accounted  for  by  incorporating  the  standard  integral- 
averaged  approach  for  the  calculation  of  the  typical  thermoelectric 
properties,  see  Sandoz-Rosado  and  Stevens  [9], 

The  resulting  system  model  consists  of  four  coupled  ordinary 
differential  equations  and  a  fifth  coupled  algebraic  Eqs.  ( 1 )— (3). 
Eq.  (1)  gives  the  thermal  governing  equations  for  two  legs,  a  and 
b,  while  Eq.  (2)  gives  the  electrical  governing  equations  for  the 
same  legs.  The  fifth  equation  (Eq.  (3))  represents  Ohm’s  Law. 


(1) 

(2) 

Mx=\)-Mx  =  \)  =  h. 

(3) 

The  complete  model  for  a  rectangular  couple  with  insulated  leg 
sides  and  two  legs  electrically  in  series,  a  and  b,  is  described  in 
terms  of  dimensionless  parameters.  The  assumed  dimensionless 
variables  of  the  model  include  a  space  coordinate  x  =  x/L  normal¬ 
ized  with  leg  length,  temperature  T  =  T /AT  normalized  with  the 
temperature  difference  driving  the  Seebeck  effect,  voltage 
<j>  =  </>/( Aa  ■  AT)  normalized  with  a  characteristic  couple  open  cir¬ 
cuit  voltage  due  to  couple  Seebeck  coefficient  Aa  =  acb  -  a„  and 
electrical  current  ?  =  IR/( Aa  ■  AT)  normalized  with  characteristic 
load  resistance  R  and  open  circuit  voltage.  Substitution  of  the 
dimensionless  parameters  into  the  governing  equation  results  in 
a  set  of  characteristic  parameters  which  govern  the  behavior  of 
all  couples  and  provide  validation  for  assumptions.  The  effect  of 
Thomson  heat  is  captured  in  the  dimensionless  parameter 
[)  =  xAa.ATL/(ARk),  which  serves  as  a  measure  of  the  accuracy  of 
neglecting  Thomson  heat  t,  in  terms  of  couple  Seebeck  Aa,  operat¬ 
ing  temperature  difference  AT,  length  L,  cross  sectional  area  A,  load 
resistance  R,  and  thermal  conductivity  k.  This  heat  is  classically 
neglected  as  part  of  the  assumption  of  constant  material  properties 
by  the  second  Kelvin  relation  t  =  T%  [1,2],  Work  in  accounting  for 
this  heat  has  been  performed  by  Sherman  et  al.  [10],  Yamashita 
[11],  and  Min  et  al.  [12],  For  f)  values  much  smaller  than  unity 
the  assumption  of  neglecting  Thompson  heat  is  valid,  but  in  the 
case  of  a  couple  with  /f  on  the  order  of  unity  this  assumption  is 
no  longer  justified.  Notice  the  term  involves  not  only  material 
properties  but  also  a  geometric  slenderness  ratio  L/A  as  well  as 
operational  parameters  AT  and  R.  The  effect  of  Joule  heating  is  cap¬ 
tured  by  the  dimensionless  parameter  y  =  Aa2  ATL2/ (A2  R2kcr),  with 
the  introduction  of  electrical  conductivity  <r.  The  dimensionless 
voltage  due  to  Seebeck  effect  is  captured  by  {  =  a/Aa,  and  the  volt¬ 
age  due  to  electrical  losses  is  found  in  X  =  L/(ARcr). 

Included  in  the  model  are  the  eight  boundary  conditions 
required  for  complete  specification  of  the  problem  statement. 
Under  a  strict  set  of  simplifying  assumptions  the  model  can  be 
seen  to  return  the  same  solution  classically  offered  as  the  analytic 
solution  of  a  thermoelectric  couple,  see  Sherman  et  al.  [10],  Lampi- 
nen  [13],  or  Bejan  [14],  These  assumptions  include:  constant  iso¬ 
tropic  material  properties,  fixed  hot  and  cold  shoe  temperatures, 
insulated  leg  sides,  steady  state,  and  one  dimensional  conduction. 
Additionally,  due  to  the  electrical  series  nature  of  the  couple  the 
electrical  current  in  legs  a  and  b  are  equal  in  magnitude  and  oppo¬ 
site  in  direction.  With  appropriate  application  of  aforementioned 
boundary  conditions  and  solution  of  the  coupled  equations  the 
classic  solution  of  the  couple  can  be  found. 

Calculation  of  the  couple’s  thermodynamic  conversion  effi¬ 
ciency  reveals  two  critical  design  parameters,  well  suited  for  opti¬ 


mization  of  device  design:  (i)  the  geometric  factor  X  =  AbLa/(AaLb) 
is  the  ratio  of  the  leg  slenderness  for  a  and  b  legs  and  (ii)  the  load 
factor  Y  =  R/(La/(Aa(Ta)  +  Lh/(Ah<jh))  is  a  ratio  of  load  resistance  to 
couple  resistance.  Details  of  the  classic  optimization  can  be  found 
in  literature  and  provide  a  means  of  calculating  the  optimized  X 
and  Y  values  to  obtain  a  couple  with  maximum  conversion  effi¬ 
ciency  [10,13-15].  In  addition  to  the  device  design  optimization, 
the  traditional  thermoelectric  figure  of  merit  can  be  extracted  from 
this  conversion  efficiency  and  used  as  a  material  design  guideline. 

Similarly  the  method  can  be  applied  to  a  cylindrical  couple  con¬ 
figured  with  radial  heat  transfer.  Legs  in  a  cylindrical  couple  are  con¬ 
figured  as  a  solid  washer  or  ring  shape  with  the  temperature  gradient 
ranging  from  the  inside  to  outside  radius.  Sets  ofp-  and  n-type  wash¬ 
ers  can  be  electrically  connected  in  series  to  complete  a  circuit 
(Fig.  1 ),  similar  to  their  rectangular  counterparts.  Cylindrical  cou¬ 
ples  may  prove  to  be  well  suited  for  the  design  of  compact  heat 
exchangers  which  require  a  radial  conduction  path.  For  instance, 
the  cylindrical  couple  would  serve  well  in  an  energy  harvesting 
application  of  a  coolant  line  passing  through  a  hot  exhaust  cham¬ 
ber.  The  outside  radius  of  the  legs  would  be  in  communication 
with  the  hot  ambient  thermal  reservoir,  and  the  cold  junction 
would  be  in  direct  thermal  contact  with  a  coolant  fluid.  Cylindrical 
couples  have  been  theoretically  and  experimentally  investigated 
by  Min  and  Rowe  [16,17],  Landecker  [18],  Lund  [19],  and  Liu 
[20],  The  previous  governing  equations  Eqs.  (1 )— (3)  can  be  substi¬ 
tuted  for  their  cylindrical  equivalents  and  solved,  see  Appendix  A. 
The  temperature  and  voltage  profiles,  now  functions  of  radius  r  and 
leg  width  w,  can  be  used  to  calculate  the  couple’s  thermodynamic 
conversion  efficiency.  Again,  two  design  parameters,  with  similar 
physical  meaning  to  their  rectangular  counterparts,  can  be 
extracted.  The  appropriate  geometric  factor  as  derived  by  this  work 
becomes  X  =  wbIn’jfL/(walnr^)  with  subscript  o  indicating 
outside  radius  and  i  inside  radius.  Likewise  the  new  cylindrical 
load  factor  derived  in  this  work  becomes  Y  =  Rj(\nlfA/ 
(2naawa)  +ln^/(27T(7f,W(,)y  Optimization  of  the  efficiency  in 
terms  of  the  cylindrical  X  and  Y  design  factors  is  identical  to  the 
rectangular  case  and  results  in  the  same  final  optimized  values, 
as  suggested  by  Min  and  Rowe  [16],  The  similarity  of  the  rectangu¬ 
lar  and  cylindrical  solutions  is  expected  and  reported  in  literature, 


Fig.  1.  Schematic  of  a  cylindrical  couple,  highlighting  the  geometric  parameters  of 
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but  the  new  cylindrical  X  and  Y  design  factors  of  this  work  are 
required  to  optimize  a  cylindrical  couple.  For  instance,  the  new 
geometric  factor  and  load  factor  allow  one  to  realize  the  maximum 
theoretical  conversion  efficiency  or  power  density  in  terms  of  cou¬ 
ple  radius  and  leg  width.  The  classic  rectangular  design  factors  do 
not  capture  geometries  useful  to  a  cylindrical  couple. 

A  useful  presentation  of  the  classic  solution,  Fig.  2,  displays  the 
behavior  of  the  couple  across  the  design  space.  The  contour  presen¬ 
tation  allows  for  a  better  physical  understanding  of  design  tolerance 
factors  over  the  simple  analytic  solution.  The  axes  of  the  plot  map 
out  a  range  of  X  and  Y  design  factors  and  the  contour  lines  mark 
the  various  levels  of  calculated  conversion  efficiency  r\.  In  addition 
to  conversion  efficiency,  lines  of  constant  power  density,  electrical  cur¬ 
rent,  or  power  may  also  be  superimposed  to  assist  in  couple  design. 
The  location  of  the  maximum  coincides  with  the  classic  analytic 
solution,  and  the  figure  provides  a  view  into  the  range  of  design 
values  which  will  return  nearly  optimized  couples.  Two  modifica¬ 
tions  were  introduced  to  this  classical  solution  with  the  new 
design  space  to  optimize  the  device  design  for  more  realistic  con¬ 
ditions.  These  modifications  include:  (i)  variable  boundary  condi¬ 
tions  and  (ii)  lateral  heat  transfer. 

Throughout  this  work  the  material  properties  employed  for  the 
models  were  selected  to  represent  NASA’s  general  purpose  heat 
source  radioisotope  thermoelectric  generator  (GPHS-RTG).  The 
Si/Ge  couples  have  an  average  n-leg  thermal  conductivity  of 
4.02  W/(mK),  electrical  conductivity  of  5.44e4  S/m,  and  Seebeck 
coefficient  of  -231  p.V/K.  The  p-leg  has  an  average  thermal  con¬ 
ductivity  of  4.24  W/(mK),  electrical  conductivity  of  4.33e4  S/m, 
and  Seebeck  coefficient  of  +208  pV/K  [21], 

2.  Variable  boundary  conditions 

In  real  life  applications  the  temperatures  of  the  hot  and  cold 
reservoirs  are  known;  the  difference  of  which  can  be  much  differ¬ 
ent  than  the  temperature  drop  across  the  legs  due  to  factors  such 
as  conductive  resistance  of  the  shoes  t/k,  interface  resistances 
^interface'  and  convection  hhc  or  radiation  resistance  between  ther¬ 
mal  reservoirs  and  the  shoes.  Therefore  in  device  design  it  is 
important  to  leave  the  thermal  boundary  conditions  of  the  legs 
free.  This  can  be  achieved  by  introduction  of  an  effective  convec¬ 
tion  coefficient. 


Conversion  Efficiency  (%) 


X 


Fig.  2.  Contour  plot  of  conversion  efficiency  for 
the  classic  assumptions. 


K?c  =  &  J  +  ^interface  +  £)  •  (4) 

This  effective  convection  coefficient  boundary  condition  can  be 
used  to  represent  an  equivalent  thermal  resistance  of  a  lumped 
combination  of  the  aforementioned  factors.  The  appropriate  ther¬ 
mal  boundary  conditions  to  apply  to  the  model  are 

-^dT^  +  Ia^fTa,m  =  fth CU  -  MO}] ,  (5) 

-ka/TjLM  +  ^Ta^)  =  h c[Ta,b(La,b)  ~  T 4,  (6) 

where  hh  and  hc  are  the  effective  convection  coefficients  of  the  hot 
(x  =  0)  and  cold  (x  =  L)  shoes  respectively  and  T„ 0|lc  are  the  ambient 
reservoir  temperatures.  The  dimensional  forms  of  the  boundary 
conditions  are  presented  for  convenience.  The  first  term  which 
accounts  for  the  conductive  heat,  and  the  last  term  which  accounts 
for  convective  heat  are  the  classic  boundary  condition  used  for  a 
non-thermoelectic  conductor,  see  Arpaci  [3],  The  Peltier  heat  is 
introduced  by  the  second  term  on  the  left  hand  side  of  the  boundary 
condition.  Constant  material  properties  and  a  negligible  Thomson 
heat  term  are  used  with  the  model.  As  discussed  previously  the 
Thomson  heat  term  is  negligible  for  constant  material  properties 
through  the  second  Kelvin  relation  [1],  Solution  using  the  full 
boundary  condition  proceeds  with  a  straightforward  uncoupling 
of  the  differential  equations,  resulting  in  a  nonlinear  expression  to 
solve  for  the  electrical  current.  The  solution  with  the  actual  bound¬ 
ary  conditions  as  described  is  solvable  but  makes  obtaining  a  closed 
form  efficiency  equation  difficult;  optimization  can  still  be  per¬ 
formed  numerically  and  will  stand  to  justify  the  simplifying 
assumption  of  neglecting  the  Peltier  heat  term  from  the  boundary. 
The  simplified  boundary  condition  leads  to  a  useful  solution  of  the 
couple,  which  allows  shoe  temperature  the  freedom  to  float  with 
the  solution  as  desired.  The  solution  introduces  a  Device  Design 
Factor,  a  dimensionless  parameter  of  range  zero  to  one,  with  the 
ideal  device  having  a  value  of  unity. 

Device  Design  Factor  :  Da,6  =  (7) 

In  the  limit  of  effective  convection  coefficients  approaching 
infinity  this  Device  Design  Factor  approaches  unity  and  the  entire 
model  reduces  to  the  classic  solution  with  the  ambient  tempera¬ 
tures  replacing  the  shoe  operating  temperatures.  For  a  real  couple 
the  Design  Factor  serves  as  an  approximation  of  the  operating  temper¬ 
ature  in  terms  of  the  ambient  temperature.  The  couple’s  realized 
temperature  difference  can  be  estimated  by  an  effective  tempera¬ 
ture,  ATEffective  =  (Da  +  D(,)AT0O/2,  an  estimation  of  the  actual  tem¬ 
perature  difference  in  terms  of  the  Design  Factor  and  ambient 
temperatures.  This  effective  temperature  can  be  used  both  to  assist 
in  couple  design  as  well  as  thermoelectric  module  integration  in  a 
system.  The  Device  Design  Factor  plays  a  significant  but  indirect 
role  in  the  couple’s  Thomson  heat  and  Joule  heating  terms  (Eq. 
(1))  by  altering  the  operating  temperatures  of  the  couple  and 
thereby  altering  the  couple’s  operating  electrical  current.  The  Fou¬ 
rier  heat  term  of  the  couple  is  also  influenced  by  the  Device  Design 
Factor  as  thermal  gradients  are  fixed  by  the  couple’s  operating 
temperatures. 

Thermodynamic  conversion  efficiency  r/  for  the  convection 
boundary  condition  can  be  derived  in  terms  of  the  Design  Factor 
D  with  the  simplifying  assumption  of  hh  =  h,  and  by  assuming  that 
the  hot  shoe  temperature  can  be  approximated  with  the  aforemen¬ 
tioned  effective  couple  temperature  (A TEff).  The  hh  =  hc  assump¬ 
tion  is  only  required  to  obtain  clean  and  meaningful  solutions  for 
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presentation;  this  assumption  is  not  required  to  hold  for  the  gen¬ 
eral  application  of  the  model. 


+  <Drt-DaS')  f1  ^  (1  D°2D *)_ 


(8) 


ZD(X,Da,D„) 


(9) 


where  ZD  is  analogous  to  the  classic  dimensionless  couple  parame¬ 
ter  Z  used  in  figure  of  merit  ZT  and  t/c=o  is  the  thermodynamic  Carnot 
efficiency  based  on  the  ambient  temperatures.  The  efficiency  can  be 
optimized  in  terms  of  the  four  design  parameters  X,  Y,  and  Da  i)  (see 
Appendix  B). 

The  efficiency  increases  with  increasing  D.  The  optimal  value  for 
D  is  unity  -  in  the  range  zero  to  unity  -  which  causes  the  conver¬ 
sion  efficiency  to  reduce  to  the  conversion  efficiency  of  the  classic 
case,  Eq.  (8).  For  a  D  of  unity  no  losses  exist  through  the  shoes  and 
the  temperature  drop  between  the  hot  and  cold  reservoirs  are  the 
same  as  the  temperature  drop  across  the  legs. 

As  a  guideline,  the  efficiency  optimized  couple  should  be 
designed  with  optimal  X  and  Y  values  in  addition  to  using  long  legs, 
having  a  high  boundary  convection  coefficient,  and  as  small  a  leg 
thermal  conductivity  as  possible.  To  exemplify  this,  the  depen¬ 
dence  of  Device  Design  Factor  is  illustrated  in  Fig.  3  as  a  function 
of  leg  length,  for  a  wide  range  of  convection  coefficients.  For  the 
clarity  of  this  demonstration  only  shoe  conduction  resistance  was 
taken  into  account,  neglecting  additional  terms  such  as  the  inter¬ 
face  resistance  and  shoe  convection.  A  shoe  convection  coefficient 
can  be  approximated  as  the  thermal  conductivity  of  shoe  material 
divided  by  the  shoe  thickness  using  Eq.  (4);  for  instance  a  1.9  mm 
thick  SiMo  shoe  would  have  roughly  13,000  W/(m2  K).  The  typical 
shoe  used  in  the  GPHS-RTG  is  1.9  mm  thick  SiMo  [21,22],  The  prac¬ 
tical  design  rule  stands  to  design  a  short  couple,  in  order  to  produce 
the  largest  power,  without  reducing  the  Design  Factor  significantly. 
For  example,  for  legs  shorter  than  10  mm  and  a  shoe  convection 
coefficient  of  50,000  W/(m2  K),  the  Design  Factor  starts  diverging 
from  unity  (Fig.  3).  In  the  case  of  1.9  mm  thick  SiMo  shoe,  with 
the  conduction  resistance  of  13,000  W/(m2  K),  legs  shorter  than 
10  mm  will  start  exhibiting  non-trivial  divergence  from  unity,  typ¬ 
ical  GPHS-RTG  legs  are  sufficiently  designed  to  be  20.3  mm  [22], 
Furthermore,  this  divergence  will  occur  at  longer  leg  lengths  as 
the  (i)  shoe  thickness  increases,  (ii)  the  conductivity  of  the  shoe 


Demonstration  of  Design  Factor  for  k=4  W/(m  K) 


material  decreases,  and  (iii)  other  factors  such  as  convection  and 
interface  resistance  are  included.  In  an  application  such  as  NASA’s 
GPHS-RTG  where  the  heat  source  is  in  direct  contact  with  the  shoe, 
the  conduction  resistance  will  dominate,  however  in  other  possible 
industrial  applications  (i.e.,  exhaust  of  automobiles)  the  convection 
will  become  the  dominant  factor  with  practical  convection  coeffi¬ 
cients  in  the  range  of  5-50W/(m2K)  [23].  In  such  applications, 
the  importance  of  this  design  parameter  becomes  even  more  sig¬ 
nificant  due  to  large  divergence  from  unity.  While  the  classic  solu¬ 
tion  predicts  the  same  conversion  efficiency  of  a  couple  of  any  length 
this  model  analytically  demonstrates  the  effect  of  length  with  the 
introduction  of  Design  Factor.  A  couple  with  Design  Factor  near 
unity  approaches  the  maximum  conversion  efficiency  as  predicted 
by  classic  theory,  while  a  reduced  factor  will  have  significant 
reductions  on  efficiency.  As  a  result  of  the  model  a  couple’s  length 
can  be  devised  in  order  to  bring  the  Design  Factor  satisfactorily 
close  to  unity,  a  reasonable  value  for  Design  Factor  may  be  selected 
to  have  a  value  such  as  0.99  so  we  can  then  calculate  the  length 
required  to  reach  this  level  as  Lg9%. 


Ld 


D(hh  +  hc)k 
(l  -  D)hhhc 


for  instance  Lg 9%  = 


99  k(hh  +  ftc) 
hhhc 


(10) 


Eqs.  (8)-(10)  provide  a  complete  description  of  the  thermoelectric 
conversion  efficiency  of  a  couple,  taking  into  account  total  thermal 
resistance  of  shoes.  The  solution  is  consistent  with  the  classic  case 
of  fixed  hot  and  cold  temperatures  in  which  convection  coefficients 
can  be  thought  of  as  being  infinite  resulting  in  Design  Factors  of 

To  this  point  the  formulation  has  neglected  the  Peltier  term  in 
the  boundary  condition,  in  order  to  obtain  the  desired  analytic  effi¬ 
ciency  equation  (Eq.  (8)).  A  solution  including  the  Peltier  term  has 
been  optimized  numerically  and  provides  a  contrast  for  the  simpli¬ 
fied  analytic  model.  A  comparison  of  the  approximate  boundary 
condition  with  the  exact  boundary  condition  is  presented  in 
Table  1.  The  table  displays  the  Xopt  and  Yopt  values  required  to 
obtain  the  optimal  conversion  efficiency  rjopt.  The  Xopt,  Yopt,  and 
ijopl  values  are  calculated  for  both  the  analytic  and  the  exact  solu¬ 
tions  over  a  range  of  boundary  convection  coefficients.  For  an  infi¬ 
nite  convection  coefficient,  Design  Factor  1.0,  the  analytic  solution 
is  identically  equal  to  the  exact  solution.  For  large  coefficients, 
greater  than  50,000  W/(m2  K),  the  difference  between  the  analytic 
model  and  the  exact  solution  are  negligible.  For  convection  coeffi¬ 
cients  below  5000  W/(m2  K),  Design  Factor  0.85,  the  difference 
between  the  analytic  model  and  the  exact  solution  becomes  appar¬ 
ent.  At  a  convection  coefficient  of  5000  W/(m2  K)  the  difference  in 
the  analytically  calculated  conversion  efficiency  (5.30%)  and  the 
exact  solution  (5.26%)  contains  less  than  a  1%  difference.  Similarly 
the  differences  between  the  analytic  and  exact  models  are  less 
than  1%  for  the  Xop,  and  Yopl  values.  Lower  convection  coefficients 
lead  to  larger  differences  between  the  analytic  model  and  the  exact 
solution.  For  a  convection  coefficient  of  500  W/(m2  K),  Design  fac¬ 
tor  0.37,  the  difference  between  the  analytic  conversion  efficiency 
(2.37%)  and  the  exact  solution  (2.28%)  is  3.9%.  The  analytic  model, 
using  the  simplified  boundary  condition,  is  found  to  be  reasonably 
close  to  the  exact  solutions.  The  simplification  is  valid  in  the  case 
osAoeA TL/(ARk)  <  1. 


3.  Lateral  heat  transfer 

The  primary  heat  transfer  in  a  thermoelectric  couple  is  down 
the  leg  length,  although  some  heat  will  travel  in  the  lateral  direc¬ 
tion  into  the  region  between  legs  or  from  the  region  between  legs 
into  the  couple.  Insulating  the  volume  between  the  legs,  or  operat¬ 
ing  in  a  vacuum  will  reduce  this  effect.  However,  it  is  of  interest  to 
investigate  the  magnitude  of  such  lateral  heat  transfer  on  the 
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Comparison  of  simplified  boundary  condition  to  numerical  solution  of  actual  condition. 


B.C.  convert 

ion  [W/(m2  K)]  Design  Factor  [-] 

Xm  (Eq.  (B.8))  [-1 

Yopt  (Eq.  (B.9))  [-] 

flopt  (Eq.  (8))  [%] 

Xopt  exact  [ 

-]  Ym  exact  [ 

-]  t]opt  exact  [-] 

oo 

1.00 

1.09 

1.22 

6.15 

1.09 

1.22 

6.15 

500,000 

0.99 

1.09 

1.22 

6.14 

1.09 

1.22 

6.14 

50,000 

0.98 

1.09 

1.22 

6.05 

1.09 

1.22 

6.05 

5000 

0.85 

1.10 

1.22 

5.30 

1.09 

1.23 

5.26 

500 

0.37 

1.11 

1.22 

2.37 

1.09 

1.27 

2.28 

thermoelectric  couple’s  operation  in  order  to  determine  critical 
design  considerations. 

The  lateral  heat  transfer  can  be  considered  using  a  one  dimen¬ 
sional  thermoelectric  conduction  equation,  derived  for  a  control 
volume  which  extends  across  the  width  of  a  leg  so  as  to  include 
heat  transfer  at  the  boundary  as  part  of  the  governing  equation. 
In  this  fashion  the  formulation  resembles  the  classic  heat  transfer 
in  a  fin  [3-5],  The  modified  thermal  governing  equations  now 
include  a  Fin  Factor,  a  dimensionless  parameter  with  the  ideal 
device  having  a  value  of  zero. 


Fin  Factor : 


f-h/l 


(H) 


where  the  leg  perimeter  P,  and  side  convection  coefficient  h  have 
now  been  introduced.  Similar  to  the  convection  boundary  condition 
solution,  the  leg  side  convection  coefficient  may  be  an  effective 
coefficient  to  account  for  any  mode  of  heat  transfer.  The  Fin  Factor 
is  found  only  in  an  additional  lateral  heat  transfer  term,  which  must 
be  included  in  the  thermal  governing  equation.  The  Fin  Factor  plays 
no  role  in  altering  the  Fourier,  Thomson,  or  Joule  heating  terms.  The 
thermal  governing  equations  now  become 


A  change  of  variables  has  been  performed  as  8(x)  =  T(x)  -  Tx 
with  T^  the  ambient  leg  side  temperature,  assumed  to  be  fixed 
along  the  leg  length  for  simplicity.  In  fact  the  ambient  temperature, 
which  the  leg  communicates  with,  may  be  a  function  of  space,  but 
the  fixed  ambient  temperature  approximation  allows  for  insightful 
investigation  without  unneeded  complexity.  The  dimensionless 
temperature  is  now  defined  as  8(x)  =  6(x)/AT.  The  remainder  of 
the  model  consists  of  Eqs.  (2)  and  (3)  in  their  presented  form  as 
the  leg  sides  are  still  assumed  to  be  electrically  insulated. 
Following  the  previous  case  the  Thomson  heat  term  will  be 
neglected  and  constant  material  properties  will  be  assumed. 
Solution  of  the  temperature  profile  is  found  to  be  hyperbolic  in 
the  space  coordinate,  in  accordance  with  a  fin  type  formulation. 
The  electrical  current  reduces  to  the  same  result  as  the  classic  case, 
as  expected. 

The  inclusion  of  an  ambient  temperature  introduces  an  inflec¬ 
tion  point,  near  the  ambient  temperature,  into  the  thermal  solu¬ 
tion,  shown  exaggerated  in  Fig.  4  by  employing  a  large 
convection  coefficient  h  =  500  W/(m2  K).  A  practical  convection 
coefficient  may  be  in  the  range  of  5-50  W/(m2  K)  [23],  The  electri¬ 
cal  current  is  fixed  to  that  calculated  in  the  classic  solution;  there¬ 
fore  the  power  output  of  the  fin  couple  is  limited  to  the  classic 
solution,  thus  for  a  system  that  is  constructed  to  have  Design  Fac¬ 
tor  of  unity.  Any  couple  which  produces  the  same  power  output  but 
draws  a  larger  heat  input  or  rejects  a  larger  heat  output  will  be  a  less 
efficient  device. 

The  four  cases  in  Fig.  4  have  hot  and  cold  end  temperatures  of 
1123  and  573  K,  respectively.  The  profiles  are  generated  with 
X  =  1  and  Y  =  1  for  convenience,  as  these  are  typical  values  used 
in  practice.  Since  the  electrical  solution  is  the  same  in  all  four  cases 
both  the  power  output  and  the  Peltier  heat  terms  are  identical. 


Typical  Temperature  Profiles  for:  X=1 ,  Y=1 ,  h=500  W/(m2K) 


Fig.  4.  Typical  temperature  profiles  for  a  range  of  ambient  temperatures,  with  X  =  1 
and  Y  =  1  for  convenience. 


Therefore  the  difference  in  conversion  efficiency  is  determined 
from  the  lateral  heat  transfer  and  the  conductive  heat  transfer 
through  the  ends  of  the  legs.  The  heat  transferred  through  the  ends 
of  the  legs  is  proportional  to  the  temperature  gradient  at  the  ends. 
For  the  extreme  case  of  an  ambient  temperature  equal  to  that  of 
the  cold  shoe,  lateral  heat  can  transfer  only  out  of  the  device  and 
a  larger  heat  input  is  drawn  for  the  same  power  output;  thus  the 
couple  will  be  less  efficient  then  the  classic  case.  From  Fig.  4  this 
case  can  be  seen  to  have  a  steeper  temperature  gradient  at  the 
hot  shoe,  compared  to  the  classic  solution,  confirming  the  larger 
heat  input.  In  the  other  extreme,  the  case  of  ambient  temperature 
equal  to  that  of  the  hot  shoe,  the  couple  will  again  produce  the 
same  power  output  and  Peltier  heat.  In  this  case  lateral  heat  is 
transferred  only  into  the  couple,  so  the  output  thermal  energy  is 
directly  proportional  to  the  thermal  gradient  at  the  cold  shoe. 
Compared  to  the  classic  solution  this  case  will  have  a  larger  heat 
output  for  the  same  power  generation,  so  the  couple  will  be  less 
efficient.  Thus  any  intermediate  temperature  can  follow  similar 
reasoning  to  be  less  efficient  then  the  classic  case.  However,  for 
intermediate  temperatures  the  lateral  heat  will  pass  into  the  leg 
along  some  portion  of  the  length,  and  out  of  the  leg  along  the 
remainder.  Therefore,  the  evaluation  of  conversion  efficiency 
becomes  more  involved;  see  Appendix  C  for  an  example.  Fig.  4 
shows  exaggerated  profiles  by  employing  a  large  convection 
coefficient. 

The  efficiency  equation  is  given  below  for  the  case  of  ambient 
temperature  equal  to  the  cold  shoe  temperature.  Two  new  geomet¬ 
ric  factors,  a  thermal  fin  geometric  factor  (G),  and  an  electrical  fin 
geometric  factor  ( H )  are  defined.  In  the  limit  of  Fin  Factor  (F) 
approaching  zero  -  the  ideal  couple  -  both  G  and  H  approach  the 
traditional  geometric  factor  X.  Although  only  Fa  is  found  in  the  effi¬ 
ciency  equation,  the  influence  of  both  leg’s  fin  factors  (F„  and  Fb) 
impact  the  efficiency  through  the  introduction  of  the  geometric 
factors  G  and  H. 
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+Y)-r,c- 


(13) 


ZF(X,C)  = 


(«b  ~  (Xq)2 

+  Gk»)  ' 


_  /PbAthtkq  tanh(Fa) 
y  PaAahaki,  tanh(F|,)  ’ 


(14) 


(15) 


jpbAhhbka  tanh(^) 

1  '  PaAahakb  tanh(^)  ’ 


(16) 


where  ZF  is  analogous  to  the  classic  dimensionless  couple  parame¬ 
ter  Z  used  in  figure  of  merit  ZT  and  t/c  is  the  thermodynamic  Carnot 
efficiency.  The  conversion  efficiency  of  the  fin  couple  simplifies  to 
the  classic  solution  when  F  approaches  zero  as  F/  tanh(F) 
approaches  unity  and  tanh(F/2) /F  approaches  1  /2.  The  design  rule 
stands  to  design  a  couple  with  a  P/A  ratio  small  enough  to  force  the 
Fin  Factor  satisfactorily  near  zero.  This  small  P/A  ratio  requires  the 
use  of  large  cross-sectional  area  legs;  for  instance  a  square  cross- 
section  leg  will  have  the  P/A  ratio  inversely  proportional  to  the  side 
length.  Therefore,  to  obtain  a  small  P/A  ratio  the  use  of  a  large  side 
length  is  required.  Using  values  typical  for  the  NASA  GPHS-RTG 
(legs  of  2.74  x  6.50  x  20.3  mm)  [22]  the  Fin  Factor  as  developed 
in  this  work  could  be  as  large  as  0.84  under  a  He  atmosphere  testing 
condition  assuming  an  h  of  5  W/(m2  K).  The  He  atmosphere  testing 
condition  is  common  for  pre-mission  tests  before  venting  to  the 
vacuum  of  space  [22],  For  non-space  applications  of  thermoelec¬ 
trics,  where  a  vacuum  is  not  available,  lateral  losses  will  always 
be  present.  While  the  classic  solution  predicts  the  same  conversion  effi¬ 
ciency  of  a  couple  of  any  P/A  ratio  this  model  analytically  demonstrates 
the  effect  of  P/A  with  the  introduction  of  Fin  Factor.  Similar  to  the  cou¬ 
ple  length  selection  of  the  previous  section,  to  keep  Device  Design 
Factor  near  unity,  a  P/A  ratio  can  be  selected  to  keep  the  Fin  Factor 
near  zero.  A  suitably  small  Fin  Factor  may  be  0.05  so  we  can  then 
calculate  the  P/A  ratio  required  to  reach  this  level  as  (P/A)5%. 


9, 


F^k 

L2hi 


for  instance 


GL 


0.05 2k 

O  ' 


(17) 


4.  Comparison  of  models 

Table  2  summarizes  a  set  of  couple  parameters  and  illustrates 
how  it  can  reduce  conversion  efficiency,  a  hazard  of  a  poorly 
designed  couple.  The  table  is  calculated  from  typical  values  for 
couples  found  in  the  GPHS-RTG,  all  cases  use  matching  couple 


material  properties  and  leg  sizes.  The  parameters  of  study  in  the 
table  are  the  end  and  side  convection  coefficients,  the  calculated 
Device  Design  Factor  and  Fin  Factor  are  also  shown.  The  effect  of 
a  range  of  Design  and  Fin  Factors  can  be  clearly  observed  on  max¬ 
imum  conversion  efficiency  and  power  density.  The  ideal  couple, 
case  1  with  Design  Factor  of  unity  and  Fin  Factor  of  zero,  gives 
the  greatest  maximum  conversion  efficiency  and  power  density. 
In  the  cases  of  Design  Factor  deviating  from  unity,  cases  2-5,  sig¬ 
nificant  reductions  in  both  maximum  conversion  efficiency  and 
power  density  can  be  observed.  Fin  Factor  is  kept  zero  for  the  study 
of  cases  2-5.  Case  2  represents  a  case  with  very  low  leg  end  ther¬ 
mal  resistance,  such  as  direct  conduction  contact,  while  case  5  rep¬ 
resents  a  condition  with  significant  leg  end  thermal  resistance, 
typical  of  a  convection  coupled  device,  combined  with  insuffi¬ 
ciently  short  leg  length.  Between  cases  2  and  5  the  conversion  effi¬ 
ciency  drops  from  6.14%  to  2.28%  and  power  density  drops  from 
17,670  to  2300  W/m2.  Additionally,  the  optimal  geometric  factor 
X  remains  unchanged,  while  the  optimal  load  factor  Y  shifts  with 
decreasing  Design  Factor.  The  reduction  in  conversion  efficiency 
in  case  5  is  the  result  of  using  the  same  couple  length  for  an  appli¬ 
cation  where  it  may  not  be  appropriate,  as  indicated  by  the  low 
Device  Design  Factor. 

The  remaining  cases  6-9  show  reduced  maximum  conversion 
efficiency  as  a  function  of  an  increasing  Fin  Factor.  Device  Design 
Factor  is  kept  to  a  value  of  one  for  cases  6-9.  As  a  result  the  max¬ 
imum  power  density  remains  constant.  The  cases  6-9  show  the 
behavior  of  the  leg  side  convection  coefficient,  with  four  orders 
of  magnitude,  a  range  which  covers  the  full  span  of  realistic 
devices.  Between  case  6  and  9  convection  coefficient  varies  from 
0.5  to  500W/(m2K)  with  a  resulting  change  in  conversion  effi¬ 
ciency  from  6.14%  to  2.70%.  The  location  of  both  the  optimal  X 
and  V  parameters  shift  as  a  function  of  the  Fin  Factor. 

The  reductions  in  conversion  efficiency  and  power  density  sug¬ 
gest  a  strong  need  to  consider  additional  design  parameters  than 
simply  the  material  figure  of  merit.  The  design  factors  introduced, 
along  with  Eqs.  (8)  and  (13),  provide  a  simple  means  of  evaluating 
couple  design  without  requiring  lengthy  work  with  numerical  sim¬ 
ulations.  Additionally,  the  parameters  of  study  introduced  in  this 
work  serve  to  help  map  out  efficient  testing  matrices  for  experi¬ 
mental  work.  For  instance  the  Device  Design  Factor  and  Fin  Factor 
incorporates  the  influence  of  several  couple  properties  and  can 
save  the  experimentalist  the  time  of  studying  the  influence  of  each 
property  individually. 

Table  2  may  serve  as  both  a  guide  for  model  verification  and  to 
assist  in  device  design.  The  change  in  conversion  efficiency  and 
power  density  as  a  function  of  Device  Design  Factor  may  be  veri¬ 
fied  experimentally  by  creating  a  set  of  couples  with  the  appropri¬ 
ate  Design  Factors  as  outlined  by  cases  2-5.  This  could  be 
accomplished  by  either  fixing  the  leg  properties  while  varying 
the  thermal  resistance  of  the  shoes  (i.e.  altering  the  shoe  thickness 
or  material),  or  alternatively  by  allowing  the  leg  length  to  vary  for 


Table  2 

Comparison  of  couple  parameters  demonstrating  the  importance  of  considering  several  factors. 


Case  #  B.C.  convection 

[W/(m2  K)] 

Leg  side  convectio 
[W/(m2  K)] 

n  Design  Factor 

H 

1-1 

Max  efficiency 
[X] 

Efficiency  opti 

mized 

Max  power  c 
[W/(m2)] 

lensity 

X„pt 

Ym 

1  oo 

0 

1.00 

0.00 

6.15 

1.09 

1.22 

17,733 

2  500,000 

0 

0.99 

0.00 

6.14 

1.09 

1.22 

17,670 

3  50,000 

0 

0.98 

0.00 

6.05 

1.09 

1.22 

17,118 

4  5000 

0 

0.86 

0.00 

5.26 

1.09 

1.23 

12,780 

5  500 

0 

0.38 

0.00 

2.28 

1.09 

1.27 

2300 

6  oo 

0.5 

1.00 

0.09 

1.09 

1.22 

17,733 

7  oo 

5 

1.00 

0.32 

6.05 

1.10 

1.21 

17,733 

8  oo 

50 

1.00 

1.00 

533 

1.20 

1.19 

17,733 

9 

500 

1.00 

3.16 

2.70 

1.59 

1.09 

17,733 
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fixed  shoe  resistance.  Through  this  study  it  should  be  verified  that 
Eq.  (10)  provides  the  minimum  leg  length  to  reduce  influence  of 
the  shoe  thermal  resistance;  thereby  providing  the  device  designer 
with  a  single  simple  guideline  for  leg  length.  Similarly  cases  6-9 
can  be  experimentally  verified  by  creating  a  set  of  couples  with 
varying  cross-sectional  area.  A  set  of  couples  could  be  created  with 
varying  leg  side  thickness,  selected  in  such  a  way  as  to  generate  a 
range  of  perimeter  to  area  ratios  which  can  satisfy  the  Fin  Factors 
outlined  in  the  table.  The  influence  of  Fin  Factor  on  conversion  effi¬ 
ciency  and  the  independence  of  Fin  Factor  on  power  should  be  ver¬ 
ified.  Additionally,  Eq.  (17)  should  be  verified  to  provide  the 
required  perimeter  to  cross-sectional  area  ratio.  The  proper  ratio 
is  needed  to  reduce  the  influence  of  lateral  heat  transfer. 


5.  Conclusion 

The  influence  of  thermal  resistance  and  lateral  heat  transfer  on 
the  behavior  of  a  thermoelectric  couple  has  been  investigated  ana¬ 
lytically.  The  analytic  approach,  as  opposed  to  numerical  simula¬ 
tions,  offers  convenient  closed  form  efficiency  equations  in  terms 
of  the  introduced  parameters  Device  Design  Factor  and  Fin  Factor. 
Using  realistic  values  of  the  Device  Design  Factor  such  as  0.38  a 
reduction  in  maximum  conversion  efficiency  from  6.15%  to  2.28% 
is  calculated.  This  factor  accounts  for  the  thermal  resistance  inher¬ 
ent  to  any  actual  thermoelectric  device.  Using  practical  values  of 
the  Fin  Factor  such  as  1.00  a  reduction  in  maximum  conversion 
efficiency  from  6.15%  to  5.33%  is  calculated.  This  factor  accounts 
for  the  lateral  heat  which  can  transfer  through  the  sides  of  a  ther¬ 
moelectric  leg.  Practical  and  simple  design  guidelines  for  the  L99% 
and  (P/A)5%  couple  parameters  have  been  generated  to  assist  in 
couple  design,  without  requiring  the  need  for  detailed  numerical 
work.  Additionally,  the  work  has  been  completed  for  both  the  com¬ 
mon  rectangular  couples  and  the  practically  interesting  cylindrical 
couples,  for  which  the  appropriate  geometric  factor  X  and  electrical 
factor  Y  have  been  derived  in  terms  of  cylindrical  parameters. 


Appendix  B.  Details  of  the  variable  boundary  condition 
solutions 


The  variable  boundary  condition  solution,  with  the  simplifying 
approximation  of  neglecting  Peltier  heat  term,  results  in  the  fol¬ 
lowing  thermal  gradient  profile 

dTa,b_  ijb  x  2kab  \  ~hh~hcAToa 

dx  ~  A2aj,Oathka/+2A2ah(7aJ<a,b  (Oa,b  V  +La,b%)  COa,b 

(B.l) 


where  coa,b  =  kab(hh  +  hc )  +  Lajlhhhc.  (B.2) 

The  variable  boundary  condition  electrical  current  solution  is 


-S2  ±  d\  -  4di  <53 
'  2Si  ' 


where  =  Cib^b  [l  -  —  ( Lbhc  +  2 kb)] - °^3- 

2Alabkh  [  0)b  2A2acak 


2 A2babkb  [ 

[l“(LA  +  2d. 


where  d2  =  -  ( .Lh  +  - 

\Ab(Jb  Aa(Ja 


R  , 


(h„hcLb 


hhcLa 


(B.4) 

(B.5) 

(B.6) 


where  <)3  =  f  — —ab - —  oc„  1  AT.*,. 

\  (Ob  COa  J 

For  small  5]  the  variable  boundary  electrical  current  approaches 


(B.7) 


The  conversion  efficiency  (Eq.  (8))  of  a  variable  boundary  couple 
is  optimized  with  the  following  parameters 
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Appendix  A.  Details  of  the  cylindrical  solutions 


The  cylindrical  thermal  equation  used  to  obtain  the  cylindrical 
form  of  the  design  parameters  X  and  Y  is  expressed  in  terms  of 
radius  r  rather  than  the  spacial  rectangular  coordinate  x.  The 
dimensional  form  of  the  governing  equation  is 


l(_k  rdT°A  ■  W  dU  4 
dr\  a,b  dr  J  2nwah  dr  4n2w2abroa,b 


(A.1) 


The  solution  of  the  cylindrical  thermal  gradient  is 


dTqjb  jj»  .  &  [(lnr,-)2-(laf#l 

dr  ^2w2abr(ra,bka,b  on2w^braaj,ka} &  I  lnrf  -  lnr0  I 


AT 

r(ln  r,  -  In  r0) ' 

The  cylindrical  electrical  current  solution  is 

(«b  -  Ota) AT 


(A.2) 


ZD(X%.Da.Dh)  = 
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(B.9) 

(BIO) 


Appendix  C.  Details  of  the  lateral  heat  solutions 

The  lateral  heat  transfer  solution  thermal  gradient  is 
=  mahChli  exp(xmab)  +  mafiC2ah  exp(-xmab),  (C.l) 


where  ma  h 


(C.2) 


where  Qab 


e  Oc  -  Oh  exp(mafiLafi) 
h  2mnh(maj>La,b) 
rexpjm^Vi,)  - 1  I 
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